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ABSTRACT 

In the present work we expand the study of time-dependent ionization previously iden- 
tified to be of pivotal importance for acoustic waves in solar magnetic flux tube simula- 
tions. We focus on longitudinal tube waves (LTW) known to be an important heating 
agent of solar magnetic regions. Our models also consider new results of wave energy 
generation as well as an updated determination of the mixing length of convection 
now identified as 1.8 scale heights in the upper solar convective layers. We present 1-D 
wave simulations for the solar chromosphere by studying tubes of different spreading 
as function of height aimed at representing tubes in environments of different magnetic 
filling factors. Multi-level radiative transfer has been applied to correctly represent the 
total chromospheric emission function. The effects of time-dependent ionization are 
significant in all models studied. They are most pronounced behind strong shocks and 
in low density regions, i.e., the middle and high chromosphere. Concerning our models 
of different tube spreading, we attained pronounced differences between the various 
types of models, which were largely initiated by different degrees of dilution of the 
wave energy flux as well as the density structure partially shaped by strong shocks, if 
existing. Models showing a quasi-steady rise of temperature with height are obtained 
via monochromatic waves akin to previous acoustic simulations. However, longitudinal 
flux tube waves are identified as insufficient to heat the solar transition region and 
corona in agreement with previous studies. 
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1 INTRODUCTION 

The outer atmospheric activity of the Sun as well 
as other late-type stars is largely determined by the 
structure and time evolution of photospheric magnetic 
fields. These fields extend into the stellar outer at- 
mosphere, where they cause nonradiative energy to 



be deposited (e.g., 


Anderson & Athav . 


19891: Linskv 


19911: Schriiverl, 2001: 


Boedan et al.. 2003: 


Roberts! 2004 


Ulmschncidcr & Musiclak, 2004: Musielak 


. I2004T). This 



energy is distributed over the chromosphere, transi- 
tion region, and corona; it is also considered pivotal 
for the heating and acceleration of solar and stel- 
lar winds. For new observations and simulations of 
the relative importance of acoustic and magnetic en- 
ergy deposition in the solar chromosphere, see, e.g., 
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iFossum fc Carlsson j2005 l). ICuntz. Rammacher fc Musielakl 
( 20071) . iKalkofenl l|2007l ) and iBello Gonzalez et alj (|2009h . 
Recent joint observational and theoretical studies aimed 
at elucidating the significance of wave heating in the 
sol ar photosphere and lowe r chromosphere were given 
bvlBeck fc Rammacher! J20ld), iBeck. Rezaei fe Puschmannl 



|2012i ), and IWedemever-Bohm et al.1 l|2012i ). The latter 
study investigates the channeling of magnetic energy 
through the solar chromosphere into the corona, thus re- 
sembling super-tornadoes under solar conditions. Notable 
reviews aimed a t solar chromospheric and co r onal heat- 
mg were given bvlNarain fc Ulmschneiderl (|l996l ). lKlimchukl 
( 2006), and lErdelvi fc Ballail (|2007l ). 



Magnetic flux tubes are a n important feat u re of 
the sola r surfa c e structure (e.g., ISpruit fc Roberts! . 1 19831 : 
ISolankil . Il996l : ISchriiver fc Zwaanl . |2000| ). Both obser- 
vational and theoretical studies showed that they ar e 
carriers of longitudinal tube waves (e.g. , iHerbold et al 



1 19851 : ISolanki fc Robertsl 1 19921 : H asan fc van Ballegooiienl . 
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2008), which give rise to considerable temperature in- 
creases as function of height as revealed by chromospheric 
spectral features. In this paper, we continue to pursue 
our previous line of research foc used on the generation 
of di fferent types of waves (e.g., iNarain fc Ulmschneiderl . 
1996), effects of the propagation and dissipation of waves 

1998 



ettects ot tne propagation and d issipation of waves 
Herbold et al.l . 1985; Fawzy. Ulmschneider fc Cunt j . 



Cuntz et all Il999h . and the emergen c e of chromo- 



spher ic emission fe.g.. [Cuntz et all Il999l ; iFawzv et al.l . 



2002; iRammacher fc Cuntz 
cord with previo us simulations 



20031 ). However, in ac- 
of acoustic waves (e.g. . 



Carlsson fc Steinl. [l9 92 . 1995; Ra mmacher fc Ulmschneiderl. 
20031 ; IRammacher fc Cuntzl f2005bl ; ICuntz et al 



2QQ7|), 



the 

described longitudinal wave simulations will also employ 
time-dependent (i.e., non-instantaneous) ionization. The 
geometrical and thermodynamic properties of the tube 
atmospheres are expected to also impact the dissipation of 
the wave energy as well as the so - called energy velocity as 
recently pointed out by IWorralll l|2012l) . Time-dependent 
ionization entails that, e.g., behind strong shocks the 
long timescales of hydrogen ionization / recombination 
initially prevent the dissipated energy to be converted into 
ionization energy, thus leading to strong temerature spikes 
as well as a variety of other dynamic phenomena. 

There exists a great motivation to revisit the dissipation 
of acoustic and magnetic waves in the solar chromosphere, 
which is the new determination of the mixing l ength n ear the 
top of the solar convective zones. IStein et all (|2009al lbh who 
provided new state-of-the-art simulations of the solar con- 
vection zone on the scale of supergranules, extending 10% 
of its depth but half of its pressure scale height, deduced a 
mass mixing length of Qm l = 1-8, thus s uperseding the pre- 
vious results of 1.5 or 2.0 i|Steffenl, Il993l ; iTrampedach et all 
which were widely used in previous solar heating com- 
putations. 

In our study we will also consider the relevance of tube 
spreading, i.e., different tube opening radii, for the energetics 
and thermodynamics of the magnetically heated solar chro- 
mospheric structure. Early results based on adiabatic longi- 
tudinal waves without the consi deration of time-dep endent 
ionization have been given by IFawzv et al.l (| 19981 ). They 
found that the tube shape is of critical importance for the 
heating of flux tubes. In fact, tubes of wide opening radii 
show little heating, whereas constant cross section tubes 
show very large heating at all heights. Previous simulations 
of longitudinal wav es for solar coronal tube structu r e hav e 
been given by, e.g., Ofman. Nakariakov fc DeForest (1 19991). 



lOfman. Nakariakov fc Sehgal ( 2000l ) . and Cuntz fc Suesa 
( 20041 ). These results further highlight the importance of 
tube spreading with respect to the shock wave amplitude 
and time-dependent and height-dependent heating rate. 
Concerning the formation of Ca II and Mg II emission, 
we will consider a two-component (magnetic and acoustic) 
model of the solar chromosphere with heating by longitudi- 
nal tube waves inside the flux tubes and heating by acoustic 
waves outside the flux tubes. 

The governing equations as well as the methods of our 
study, including the computation of the initial flux tube 
models, are discussed in Sect. 2. In Sect. 3, we describe the 
results of our longtitudinal wave simulations for solar flux 
tube models with different tube spreadings; the latter cor- 
respond to different magnetic filling factors. The focus of 



these studies concerns the effects of time-dependent ioniza- 
tion. Our summary and conclusions are given in Sect. 4. 



2 METHODS 

2.1 Theoretical approach 

In the following, we summarize some of the key equations 
with focus on longitudinal MHD tube waves propagating 
along the vertically dir ected magnetic flux t ubes. Following 
previous work by, e.g.. iHerbold et ail (JT985) , the equations 
to be solved in the Euler frame for continuous 1-D flows are 
given as 



dpA dpuA 
dt dx 
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rad 
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(3) 
(4) 
(5) 



Here all variables have their usual meaning. Note that p 
and p e are the gas pressures inside and outside the tube, 
respectively, B is the magnetic field strength, and e = p e : p 
denotes the pressure ratio between outside and inside of the 
tube. Furthermore, the cross section A is time-dependen t 
owing to the distensibility of the tube (e.g.. iLighthilll . Il978l ) . 
Moreover, Eqs. (|4]) and ([5]) represent the conservation of the 
magnetic flux <E> and the horizontal pressure balance between 
inside and outside the tube. All quantities are both functions 
of height x and time t; strictly speaking, this also applies to 
p e where we howe ver assume p e = p e (x) in accord to pr e- 



vious studies (e.g.. lHerbold et all H985I ; ICuntz et ail . 1 19991 ). 
as we do not consider the effects of the time-dependency of 
the outside medium toward the tube. 

Similar to t he case of acoustic waves (e.g. , 
Ulmschneider et all 1 19771 ; ICuntz fc Ulmschneiderl. 

of equations can be transformed 
form as pursued in detail by 



19881 ). the system 
into its characteristic 



IRammacher fc Ulmschneiderl (|2003l ). A key equation con- 
cerns the behaviour of the thermodynamic variables along 
the C° characteristics with the relevant compatibility 
relation given as 



dS = ^ 
dt 



dt 



which can also be written as 



./-if: 

dp 



dp 



/dT\ dS_ 
\dSj v dt 



dt = 



(6) 



(7) 



Note that, e.g., (dT/dp) s and (dT/dS) p consti- 
tute well-known therm odynamic relationships (e.g., 
iMihalas fc Mihalad . I1984T ); they need to be solved numeri- 
cally if deviations from the thermodynamic equilibrium, as 
in case of our study, are considered. 

The rate equations for the general case of flows have 
been discussed in detail by, e.g., IMihalas fc Mihalasl (|l984l ) 
and have been implemented into actual OHD/MHD codes 
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by ICarlsson fc Steinl <|2002l) . | Rammacher fc Ulmschneiderl 
(2003). among others. Assuming a simplified hydrogen atom 
of TV = 3 levels, namely two bound levels plus the continuum 
level, while considering radiative processes R and collisional 
processes C, the population density of level i can be obtained 
through solving 



drii driiii 
dt dx 



N 

E 

= 1 A j^i 



rijPji 



where Pit denotes the rate of transitions (per cm 3 and s) 
taking place from level I to level k and Pik = Rik + Gih- 
This equation represents the conservation equation for the 
particle number density m; however, it also takes into ac- 
count the generation and destruction of m caused by transi- 
tions from and to other levelfl The LHS of Eq. © would be 
zero if the particles of species i were conserved; in this case, 
Eq. ([8]) would represent the standard continuity equation for 
particles of species i. 

If time-dependent (i.e., noninstantaneous) ionization 
processes are considered, the change in the population den- 
sity for the continuum level is given as 



dn c 
~~dt 



N-l N-l 

^ m(Ric + C ic ) - n c ^2(R cl + Coi) - n c W (9) 



(Rammacher & Ulmschneiderl . 120031 ): the latter also requires 
the time-dependent solution of the advection term given as 



W 



8T 

dp 



dS 
dt 



1 dp 

pc| dt 



(10) 



where cs denotes the adiabatic sound speed given as 



f's 



-, 1/2 



or 
dp 



1 -1/2 



T + (ar) p v _, /s 

The number density of the hydrogen continuum pop- 
ulation is closely related to the electron density n e , 
which is of critical importance for the manifestation 
of the chromospheric rad i ative e mission functi on (e.g., 
IVernazza. Avrett fc Loeserl . Il98ll ; lAvrettl Il985l ). Other 
equations concern the conservation of the total number of 
hydrogen particles n to t, given as 



(11) 



JV-l 

1 + Z e i 



1=1 



Tli + 



/ P 



and the conservation of electrons, given as 

N-l 



He 



+ 



^2 Hi , 



= ntot , (12) 



(13) 



where X a \ and Z a \ denote the element abundances of hy- 
drogen and the metals, respectively. Note that p is the gas 



1 Equation l[8|l is in li n e w ith the earlier wo r k by 
I Rammacher fc Ulmschneiderl ll2003f ) . In ICarlsson fc Steinl (|2002} l 
based on a 6-level atom it was shown that hydrogen ionization 
goes with collisional excitation to n = 2, then with radiative 
ionization to the continuum, a pattern broadly consistent with 
the present study. On the other hand, ICarlsson fc Stem] l|2002l 'l 
found that recombination transitions occur by going from the 
continuum level to excited levels with n > 2, which is not 
possible in the 3-level atom selected. 



pressure, T is the temperature, and k is Boltzmann's con- 
stant. 



2.2 Comments on the (magneto-)hydrodynamic 
computer code 

The treatment of acoustic and longitudinal flux 
tube waves concerning solar and stellar chromo- 
spheric models is pursue d based on the code by 
Rammacher fc Ulmschneider (l2003h and subsequent 



augmentations bv lRammacher fc Cuntzl l|2005bl ). This code 
is suitable for the simulation of one-dimensional wave 
propagation and dissipation and, moreover, also allows 
the treatment of time-dependent (i.e., non-instantaneous) 
ionization processes for hydrogen and other elements 
by obtaining solutions of the time-dependent st a tistica l 
rate equations (see also Carls son fc Steinl 1 19921 . Il995l ). 
Additional features of the code encompass the evaluation 
of the radiative losses (and gains) by also taking into 
account departures from local thermodynamic equilibrium 
(NLTE). Shocks are treated as discontinuities based on 
adequate solutions of the Rankine-Hugoniot re l ations ; 
see, e.g., previous work by iNieuwenhuiizen et al.l (|l993l ) 
for a quasi-analytic method for the computation of the 
thermodynamic relationships across shocks in the presence 

of radiative and ioni zation processes. 

The method of iRammacher fc Ulmschneiderl l|2003l ) is 
suitable for the treatment of both monochromatic and spec- 
tral waves. In the framework of our models, time-dependent 
ionization is considered with respect to hydrogen as well 
as magnesium and calcium. Regarding the chromospheric 
emission losses, we provide a detailed treatment of the 
Ca II K and Mg II k lines; the total chromospheric emis- 
sion losses a r e obt ained via appropriate scaling; see, e.g., 
ICuntz et all (l999) for previous applications of this ap- 
proach. The scaling factors have been determined by in- 
specting representative solutions of the wave models and 
by applying multi-level radiative transf er computations with 
MULTI (see IRammacher et all 120051 ). Radiative transfer 
in the con tinuum, notable H~, is pursued f o llowin g the 
method of ISchmitz. Ulmschneider fc Kalkofenl |l985). An- 
other important aspect of our calculations is the imple- 
mentation of boundary conditions. At the top boundary, 
a transmitting boundary condition is used, whereas at the 
bottom boundary, the injection of the magnetic or acous- 
tic wave energy flux is simulated via a piston-type bound- 
ary condition; the latter is suitable for monochromatic 
waves as well as acoustic and LTW frequency spectra; see 
IRammacher fc Ulmschne ider (2003]) and references therein 
for details. 



2.3 Wave energy fluxes and wave energy spectra 

The wave energy fluxes and wave energy spectra are cal - 
culat ed following the approach by Musielak et al.l (|l994l ) 
and [ulmschiieider. Theurer fc Musielakl ( 19961 ); see also 
iFawzv fc Cuntzl (|201ll ) for updated results pertaining to a 



grid of models for a set of main-sequence stars including the 
Sun. These types of models incorporate a detailed descrip- 
tion of the spatial and temporal spectrum of the turbulent 
flow to obtain adequate results for the frequency integrated 
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Table 1. Acoustic and LTW energy generation. 
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Figure 1. Depiction of the generated flux for acoustic waves 
(solid) and longitudinal flux tube waves (dashed) as function 
of frequency for a mixing length of oml = 1-8. For the mag- 
netic waves we consider the model with a bottom magnetic field 
strength of 1700 Gauss. 



acoustic energy fluxes along with the wave frequency spec- 
tra; furthermore, they utilize an extended Kolmogorov spec- 
trum with a modified Gaussian frequency factor. Both the 
acoustic and magnetic wave energy fluxes depend on the 
mixing-length para meter oml, for which «ml = 1.8 is used 
(|Stein et all 12009 d lbh. 

According to the updated value of ckml , the initial wave 
energy flux for the acoustic model is given as Fm = 1-09 x 
10 8 erg cm" 2 s" 1 (|Ulmschneider et all GUjj ); see Table 1. 
For the wave energy flux of the longitudin al flux tube wave 
we u se Fm = 2.80 x 10 8 erg cm" 2 s" 1 |Fawzv fc Cunt j . 
l201lfi : this value assumes a pressure ratio between outside 
and inside of the tube of e = 3 (see Eq.[S]). It is based on the 
assumption of an inside magnetic field strength of 1700 G, 
noting that in our model the equipartition magnetic field 
strength is given as B cq — ^/Wirpl = 2082 G with B/B eq ~ 
0.82. The adopted inside tube magnetic field strength has 
been taken at an optical depth of T5000 = 1-3, corresponding 
to a mass column density of 2.1 g cm -2 . It is found that 
the longitudinal wave energy flux is increased relative to the 
acoustic energy flux by a factor of 2 to 3.5, depending on the 
mixing- length parameter oml, owing to decisive differences 
in the turbulent velocity vt, which is considerably enhance d 
in the magnetic case (e.g. . lUlmschneider fc Musielakl .ll998). 

For previous elabor ations on the appropriate value 
for Fm and B e q see, e.g., lUlmschneider. Musielak fc Fawzvl 
|200ll ). Detailed studies of magnetic field s trengths for so - 
lar flux tubes have been given by, e.g., IStenflol (|l978l ) 



and ISolankil 1 19931 ), allowing us to motivate our choice 
of B Bq . Information on the wave frequency spectra for 
longitudinal flux tube w aves was presented by, e.g., 
lUlmschneider fc Musielakl (| 19981 ) and lUlm schneide r et al.l 
|200ll ); see Fig. 1 for the spectral wave energy distribution 
serving as basis for our modelfl These authors found that 
the spectral energy distribution for both the acoustic and 
longitudinal frequency spectrum has a maximum close to a 
wave period of 60 s. As our study also includes monochro- 
matic waves, this value will be used in the following for both 
the acoustic and longitudinal wave study. 

There is also considerable previous work on the most 
appropriate value for the initial wave energy flux of longi- 
tudinal tube waves; the latter is also modestly affected by 
the solar photospheric opacities owing to their influence on 
the construction of the solar tube models (i.e., attainment of 
radiative equili brium, position of optica l depth rsooo). Pre- 
vious work bv lUlmschneider et al. |2001 ) adopted the opac- 
ity tab les compiled bv lBohnl (|1984[ ) and lUlmschneider et al.l 
(1996), whereas more recent simulations (|Fawzvl . |2010| ) use 
the opacity table given by R . L. Kurucz and collaborators 
(see lCastelli fc Kurucj|2004l for details). 

This latter approach is adopted in the present paper. 
In principle, the opacity table by R. L. Kurucz and collab- 
orators yields noticeably lower initial wave energy fluxes, 
which for «ml between 1.8 and 2.0 are reduced by typi- 
cally 30% compare d to the models based on the opacity 
table considered bv lUlmschneider et all (|200lh . A relatively 
high initial wave energy flux base d on om l = 2.0 was also 
adopted by iRammacher fc Cu ntz ( 2005a|) , which is a fur- 
ther reason for revisiting the propagation and dissipation of 
longitudinal tube waves in the Sun by considering an ad- 
vanced treatment of the hydrodynamic and thermodynamic 
features as well as a realistic initial wave energy flux. Fu- 
ture investigations of stellar convection and photospheres 
may also be based on the massively para llel recently devel- 
oped code Bifrost (|Gudiksen et al.l . l201ll ). which is designed 
to simulate stellar atmospheres from the convection zone to 
the corona. 



2 The longitudinal wave energy spectrum as obtained through 
the study of nonlinear time-dependent responses of theoret- 
ical solar flux tubes to exter nal pressure fluctuations (e.g., 
Ulmschncidcr & Musielak, 1998) is found to be somewhat bumpy, 
which is a consequence of the occurrence of large-amplitude per- 
turbations (i.e., spiky waves) that are still apparent after temporal 
averaging has been applied. 
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Figure 2. Spread of magnetic flux tubes for the different mag- 
netic filling factors. 



Table 2. Magnetic tube models. 







r T 


B(rr) 


(%) 


(km) 


(km) 


(G) 


0.1 


55 


150 


229 


0.01 


55 


700 


10.5 


0.001 


55 


1500 


2.3 



of 700 and 1500 km is a factor of 4.6. The flux tubes repre- 
sent regions with magnetic filling factors / = r\ jr\ equal to 
0.1%, 0.01% and 0.001%, respectively, with r-g, and tt denot- 
ing the base radius and top radius of the tube, respectively. 
The effect of flux tube spreading results in a progressively 
lower inside magnetic field strength with a decreasing mag- 
netic filling factor /, entailing B(tt) = 2.3 G for / = 0.001 
compared to 229 G for / = 0.1 (see Table 2). 



1.5 
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50 100 150 200 250 300 350 400 450 500 
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Figure 3. Velocity at the bottom boundary of the flux tube mod- 
els in response to the longitudinal wave energy spectrum. 

2.4 Flux tube models 

For our models we consider three different flux tube geome- 
tries characterized by d ifferent amounts of spreading as func - 
tion of height; see , e.g..lRammacher fc Ulmschneiderl (|l989l ) 
and iFawzv et all (|l998l ) for previous work on the effects of 
tube spreading on the deposition of non-radiative energy 
and the formation of chromospheric emission with focus on 
the Sun. The bottom radii of the tube models are chosen 
as 55 km corresponding to half of the photospheric pres- 
sure scale-height; note that the base press ure of the tubes 
is clo sely aligned to the VAL-C model l|Vernazza et all 
ll98lT >. We also assume a magnetic field strength at the bot- 
tom of the tubes of Bp = 1700 G informed by previous 
studies (|Fawzv et all 1 19981 ). We again consider wine-glass 
shaped tubes with top opening radii of 150 km, 700 km 
and 1500 km, respectively (see Table 2). For these tubes, 
exponential spreading is assumed at relatively low heights, 
followed by linear spreading. Full spreading is attained at 
heights of 571 km, 1126 km and 1639 km, respectively; see 
Fig. 2 for a depiction of the adopted tube models. 

The tube with a top opening radius of 700 km has a 
top opening area approximately twenty times larger than 
the tube with 150 km. Moreover, the difference in the open- 
ing areas between the wine-glass tubes with opening radii 



3 RESULTS AND DISCUSSION 

3.1 Detailed time-dependent model simulations 

By adopting the wave energy flux and wave frequency spec- 
trum for «ml = 1.8, we pursued detailed simulations for 
longitudinal tube waves for three different flux tubes (see 
Sect. 2.4); in the following, we show calculations pertain- 
ing to the magnetic filling factors of / = 0.1% and 0.001%. 
The stochstic velocity function at the bottom of the tubes 
owing to the application of the longitudinal wave frequency 
spectrum is depicted in Fig. 3. Time-dependent ionization 
has been considered for hydrogen as well as for Ca II and 
Mg II. To demonstrate the structure of our dynamic models 
in conjunction with the impact of different tube spreadings, 
we depict snapshots at identical elapsed times, which are: 
400, 800, and 1400 s (see Fig. 4). 

The general behaviour of the depicted waves in these 
models is characterized by temporal increases of the wave 
amplitudes resulting in the formation of shocks, typically 
occurring near 700 and 800 km for / = 0.1% and 0.001%, 
respectively. The tube geometry determines the general po- 
sition of the chromospheric temperature rise. It is found that 
the larger the tube spreading, the greater the height of the 
temperature rise. In models of negligible tube spreading, the 
shock amplitudes tend to stay constant or increase as func- 
tion of height. However, this behaviour is counteracted by 
the impact of tube spreading: the larger the spreading of 
the tube, the greater the area over which the wave energy is 
distributed (i.e., geometrical dilution), resulting in smaller 
wave amplitudes (see also Sect. 3.3). 

Both models are shaped by strong interaction of shocks, 
which leads to increased shock strengths M s as function of 
height. The shock strength M s is given as 



CTl 

where (7 s h denotes the shock speed, and ui and cti denote 
the gas speed and tube speed in front of the shock, respec- 
tively. 

This behaviour is most pronounced for relatively narrow 
tubes; for example, we find that the main shock at t — 1400 s 
in the model of / = 0.1% has a strength of 6.62 compared to 
1.02 in the model of / — 0.001%. This type of behaviour is 
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Figure 4. Two sequences of snapshots of longitudinal wave computations for magnetic flux tubes with filling factors of / = 0.1% (left 
column) and / = 0.001% (right column) based on time-dependent ionization and with consideration of frequency spectra. These filling 
factors correspond to tube opening radii of 150 km and 1500 km, respectively. We depict the following variables: temperature (solid 
lines), gas speed (dashed lines), density (short-dashed lines), sound speed (long-dashed lines), and gas pressure (dashed-dotted lines). 
The snapshots are taken at times of 400, 800, and 1400 s (top to bottom). 



also mirrored by the jump in temperature across the shocks. 
The corresponding post-shock temperatures are typically 
considerably higher at relatively large atmospheric heights 
and in models of relatively narrow tubes, noting that nar- 
row tubes imply higher shock strengths anyhow and, more- 
over, entail more intense shock merging. For / = 0.1% at an 
elapsed time of 1400 s (see Fig. 4), the attained post-shock 
temperatures range between 4650 and 16,170 K. In models 
of smaller magnetic filling factors, the post-shock tempera- 
tures are typically considerably lower; additionally, the vari- 
ations in post-shock temperatures for given models simula- 
tions and time steps is also reduced. Furthermore, in models 
of time-dependent ionization, the post-shock temperatures 
are also considerably higher than in equivalent models with- 
out time-dependent ionizatio n, a behaviour consistent with 
previous acoustic models (e.g.JCarlsson fc Steinl ["l992. 1995; 
iRammacher fc Ulmschneiderl 120031 ); see also Sect. 3.2. 

Additionally, the occurrence of stronger shocks in mod- 
els of / = 0.1% compared to / = 0.001% results in higher 



levels of momentum transfer, which typically lead to smaller 
densities in the top regions of those tubes as well as charac- 
teristic behaviours for the time-averaged temperatures (see 
Sect. 3.3). Although relatively large gas speeds may be at- 
tained by this process, the outflow speeds typically remain 
subsonic. Additionally, due to energy dissipation and radia- 
tive energy losses, the wave energy fluxes considerably de- 
creases with atmospheric height (see Table 3 and 4). How- 
ever, comparing tubes of / = 0.1% and 0.001%, the decrease 
in the wave energy flux with height regarding the tube with 
/ = 0.001% is significantly less than expected from the de- 
gree of geometrical dilution. This phenomenon is caused by 
the relatively low density in the tube with / = 0.1% owing 
to the action of strong shocks (see Sect. 3.3); for previous 
result s on adiabatic LTW simulations see, e.g.. iFawzv et all 
(1998). At a height of 2000 km, the wave energy fluxes are 
given as 5.10 x 10 3 and 1.93 x 10 3 ergs cm -2 s _1 for / = 0.1% 
and 0.001%, respectively (see Table 4). For the intermedi- 
ate models of / = 0.01% (not shown), the wave energy 
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Figure 5. Behaviour of the temperature (left column) and hydrogen ionization degree (right column) of propagating longitudinal flux 
tube waves concerning tube models with different spreadings, which correspond to filling factors of / = 0.1%, 0.01%, and 0.001%. We 
depict monochromatic wave models (top) and models with frequency spectra (middle) without time-dependent ionization as well as 
models with frequency spectra including time-dependent ionization (bottom). Note the vast differences in the temperature amplitudes. 
The very large temperature spike in the model with / = 0.1% (bottom) reaches 16,170 K. 



flux at 2000 km is given as 4.17 x 10 3 ergs cm -2 s -1 . At 
low chromosphere heights, the heating and energy dissipa- 
tion rates are identified as almost independent of the mag- 
netic filling factor cons istent with observational constraints 
l|Bruls fc Solankil . ll993l ). 



3.2 Effects of flux tube spreading and 
time-dependent ionization 

To obtain further insight into the effect of tube spreading 
on the atmospheric shock strengths as well as into other fea- 
tures owing to time-dependent ionization we computed a set 
of detailed models. They include monochromatic wave mod- 
els without time-dependent ionization as well as models with 
frequency spectra with and without time-dependent ioniza- 



tion. For each set of models, we considered tube spreadings 
corresponding to magnetic filling factors of / = 0.1%, 0.01%, 
and 0.001%, respectively, thus allowing us to gauge the im- 
pact of tube spreading on the temperature structure and 
shock strengths. 

The first set of models is based on monochromatic 
waves. For all three tube spreadings we chose an elapsed time 
of 2500 s, thus ensuring appropriate comparisons between 
the models. Hence, almost 30 shocks have been inserted 
into the atmosphere; therefore, all three tube atmospheres 
have reached a dynamic steady-state, i.e., the interaction be- 
tween the shocks (i.e., ongoing shock merging) has subsided. 
Comparions between models of different tube spreading (see 
Fig. 5; top panel) indicate that the models with narrow 
tube spreading (i.e., high magnetic filling factor) are charac- 
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Table 3. LTW energy flux: monochromatic waves." 



Height 




nTDI 




TDI 






/ = 


= 0.1 


/ = 0.001 


/ = 


= 0.1 


/ = 


0.001 


(km) 
















600 


1.43 


xlO 7 


7.48 xlO 6 


1.08 


xlO 7 


1.56 


xlO 6 


900 


2.21 


xlO 6 


4.15 xlO 5 


1.97 


xlO 6 


1.01 


xlO 5 


1200 


2.62 


xlO 5 


1.01 xlO 4 


2.56 


xlO 5 


2.79 


xlO 3 


2000 


1.02 


xlO 4 


3.37 xlO 3 


4.87 


xlO 3 


2.67 


xlO 2 



a The LTW energy fluxes are given in ergs cm 2 s . 



Table 4. LTW energy flux: spectral waves. a 



Height 




nTDI 




TDI 






/ = 


= 0.1 


/ = 


0.001 


/ = 0.1 


/ = 


0.001 


(km) 
















600 


1.09 


xlO 7 


3.24 


xlO 6 


1.53 xlO 7 


3.49 


xlO 6 


900 


2.75 


xlO 6 


4.90 


xlO 5 


3.33 xlO 6 


5.49 


xlO 5 


1200 


4.30 


xlO 5 


3.31 


xlO 4 


5.07 xlO 5 


2.58 


xlO 4 


2000 


4.34 


xlO 4 


5.41 


xlO 3 


5.10 xlO 3 


1.93 


xlO 3 



a The LTW energy fluxes are given in ergs cm 2 s . 



terized by relatively high shock strengths, whereas models 
with wide tube spreading (i.e., low magnetic filling factor) 
are shaped by relatively small shock strengths. Specifically, 
for regions beyond 1200 km, the typical shock strength for 
/ = 0.1% varies between M s — 1.85 and 2.77, whereas for 
/ = 0.001%, it is only M s = 1.60. For / = 0.01%, inter- 
mediate shock strength values are found, as expected. In 
models of monochromatic waves the shocks have attained 
limiting shock strength, which is found to depend on the 
height-dependent behaviour of the tube sprea ding, a s prev i- 
ously pointed out through analytical means l|Cuntj.[200i ). 
Previous simulations for solar flux tubes based on adiabatic 
LTW waves without time- dependent ionizatio n also show a 
similar kind of behaviour (|Fawzv et al. 1. 1 19981 ). Due to the 
omission of detailed time-dependent ionization, it is found 
that the degree of hydrogen ionization at 2000 km is about 
95% for / = 0.1% and 0.01%. For / = 0.001%, a lower value 
for the hydrogen ionization degree is found; in this case, a 
smooth and steady increase of the hydrogen ionization de- 
gree occurs between 1200 and 2000 km. 

The second and third sets of models are based on LTW 
frequency spectra. Figure 5 (middle and bottom panel) 
shows snapshots without and with the consideration of time- 
dependent ionization, respectively. For both types of mod- 
els, we again consider magnetic filling factors of / = 0.1%, 
0.01%, and 0.001%. The snapshots displayed for these sets 
of models are taken at 800 and 1400 s, respectively; fur- 
thermore, the elapsed time of simulation is about 2000 s. 
The main shocks in Fig. 5 (middle panel), attained through 
the process of shock merging, have strengths of 3.19, 3.18, 
and 2.11, respectively. The main shock for / = 0.1% has a 
post-shock temperature of 10,070 K; the hydrogen ionization 
degree in front of the shock is 78%, and shortly behind the 
shock it is close to zero. By comparison, the post-shock tem- 
peratures of the main shocks pertaining to / = 0.01% and 
0.001% are 9250 K and 6470 K, respectively, which clearly 



indicates the impact of tube spreading. Clearly, the change 
of structure due to different magnetic filling factors is a di- 
rect consequence of the dilution of the wave energy flux (see 
Table 4), including associated magnetohydrodynamic, ther- 
modynamic and radiative phenomena, particularly those as- 
sociated with the formation of shocks of different strengths. 
At heights of 1200 km, the wave energy flux is reduced by 
factors of 1.5 x 10~ 3 , 2.9 x 10" 4 , and 1.9 x 10" 5 in tube 
models with / = 0.1%, 0.01%, and 0.001%, respectively. 

Figure 5 (bottom panel) reflects a similar setting but 
now both time-dependent ionization and LTW frequency 
spectra are taken into account. The main difference to the 
kind of simulations previously discussed is an enhanced ten- 
dency of building up very strong shocks due to shock merg- 
ing, particularly for the model with a magnetic filling factor 
of 0.1%. In this case, a very large temperature spike with 
a post-shock temperature of 16,170 K is found to occur; 
the corresponding shock strength is given as M s — 6.62. In 
contract, however, at relatively low heights, a considerable 
number of small shocks is encountered. Moreover, at heights 
between 600 and 1000 km, shock strengths are typically be- 
tween M a — 1.2 and 1.5 in models of / = 0.1%. The shock 
strengths are even lower in models with 0.001%; in the lat- 
ter cases, the shock formation is postponed to larger heights. 
The merging of shocks leads to the build-up of very strong 
shocks, a process that is significantly amplified in models of 
time-dependent ionization, owing to the fact that virtually 
no energy is used immediatelly behind the shocks to ion- 
ize hydrogen (or any other species). Therefore, an increased 
amount of energy is available to further increase the post- 
shock temperatures. 

The dilution of the wave energy flux is also highly sig- 
nificant in these types of models depending on the degree 
of geometrical spreading. At heights of 1200 km, the ini- 
tial wave energy flux is reduced by factors of 1.8 x 10~ 3 , 
3.9 x 10" 4 , and 9.2 x 10" 5 , in models with / = 0.1%, 0.01%, 
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Figure 6. Build-up of mean (i.e., time-averaged) temperatures 
for models with and without time-dependent ionization (bottom 
and top figure, respectively). Note the influence of the different 
tube spreadings, denoted as magnetic filling factors / = 0.f% 
and 0.0f%; see Table 2. The time of averaging for the various 
models is close to 60 s. For the simulations with time-dependent 
ionization the averaging starts at 5f3 s (/ = 0.1%) and 821 s (/ = 
0.01%) and for the simulations without time-dependent ionization 
it starts at 607 s (/ = 0.1%) and 612 s (/ = 0.01%). 



and 0.001%, respectively. However, an assessment of the ab- 
solute amount of wave energy flux indicates (see Table 3 
and 4) that it is less diluted than expected from the degree 
of geometrical dilution, which obviously is a factor of 100 
between / = 0.1% and 0.001%. The reason for this type of 
behaviour is that models with / = 0.1% have a lower overall 
density; furthermore, in those models there is a higher loss of 
energy behind the shocks due to enhanced radiative energy 
losses and the initiation of episodic outflow. For the same 
reason, the reduction of the wave energy flux with height is 
also more drastic in models with time-dependent ionization 
compared to models without consideration of this process 
(TDI versus nTDI). 

For the assortments of small shocks, it is noteworthy 
that different shock strengths correspond to different shock 
speeds; therefore, the shocks in models with LTW frequency 
spectra are poised to overtake one another, thus invoking 
significant shock merging as well as the build up of very 
strong shocks. This resu lt has already been obtained in 1- 
D acoustic models (e.g.. ICuntj . 119871 ). and it is evident in 
1-D LTW models with and without the consideration of 
time-dependent ionization, but it is drastically enhanced if 



time-dependent ionization of hydrogen is considered. The 
occurrence of very strong shocks is a typical feature of the 
combined influence of frequency spectra and time-dependent 
ionization i nherent in 1-D (magneto -) hydrodynamic mod- 
els; see, e.g., ICarlsson fc Steinl |l992t) for previous results on 
acoustic models. Figure 5 (middle and bottom panels) also 
depicts the behaviour of the hydrogen ionization degree in 
the various models in front and behind the shocks. Clearly, 
the hydrogen ionization degrees are progressively lower in 
models of smaller values of /; however, in models without 
time-dependent ionization of hydrogen it drops to almost 
zero in the post-shock regions (see Fig. 5, middle panel), 
whereas in models with time-dependent ionization of hydro- 
gen it is essentially maintained (see Fig. 5, bottom panel). 



3.3 Behaviour of time-averaged temperatures and 
ionization degrees 

Additional insight into the overall dynamic structure of flux 
tubes subjected to the propagation and dissipation of lon- 
gitudinal tube waves can be obtained by assessing the be- 
haviour of time-averaged temperatures. Again, we consid- 
ered both monochromatic and spectral waves, while pursu- 
ing detailed comparisons for models with and without time- 
dependent ionization. Concerning monochromatic models, 
we also pursued computations for acoustic waves aimed 
at representing atmospheric structure exterior to the flux 
tubes; in this regard, plane-parallel geometry was assumed. 
In all types of models, the wave energy generation was de- 
rived from a mixing length of qml = 1.8 resulting in an ini- 
tial wave energy flux of 1.09 x 10 s and 2.80 x 10 s erg cm -2 s _1 
for the acoustic and magnetic models, respectively (see Ta- 
ble 1). For the magnetic flux tube simulations, we again 
considered tube geometry given by magnetic filling factors 
of / = 0.1%, 0.01%, and 0.001%. Examples for the steady 
build-up of time-averaged temperatures owing to the imple- 
mented averaging procedure are conveyed in Fig. 6. 

Detailed information on the behaviour of the time- 
averaged temperatures is given in Fig. 7, which depicts a 
total of 14 models. Let us first focus on the evaluation of 
dynamic magnetic flux tube models based on monochro- 
matic waves. Typically, we awaited the insertion of 20 to 
50 wave periods into the atmosphere prior to the evalua- 
tion of mean (i.e., time-averaged) quantities to ensure that 
the switch-on behaviour of the tube atmospheres has sub- 
sided and dynamic equilibria have been reached. The time- 
averaged temperatures of the various tube models as well as 
those of the acoustically heated models are relatively sim- 
ilar at heights below 800 km, although they are found to 
be highest in the model with the smallest spreading (i.e., 
/ = 0.1%), as expected. This latter result can be under- 
stood based on the dilution of the wave energy flux, which 
also affects the height of shock formation. The latter is given 
as 400, 550, and 850 km for the models with 0.1%, 0.01%, 
and 0.001%, respectively. For the acoustic model, the height 
of shock formation is given as 600 km. Only minor differ- 
ences are obtained between models with and without time- 
dependent ionization, noting that due to the relatively high 
densities, low temperatures and small shock strengths (if 
any), the effects due to time-dependent ionization are less 
pronounced. 

In the upper parts of the tube atmospheres, the tem- 
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peratures in the model of the widest tube opening radius 
are expected to be lowest due to the largest dilution of the 
wave energy flux. The difference in dilution is expected to 
be a factor of 100 between / = 0.1% and 0.001%; how- 
ever, it is significantly lower, i.e., a factor between 3 and 
26 for models with or without time-dependent ionization 
and with or without the consideration of LTW frequency 
spectra (see Table 3 and 4). For some portions of the tube 
atmospheres, it is even found that "the larger the tube open- 
ing radius, the higher the mean temperatures" , a result that 
is highly counter-intuitive. The reason for this behaviour is 
that the dilution of the wave energy flux does not exactly 
follow geometrical scaling. Other relevant processes involve 
the density structure within the flux tubes as well as signif- 
icant dynamical and radiative processes, especially pertain- 
ing to strong shocks, including significant radiative energy 
losses. Comparing models with and without time-dependent 
ionization reveals that time-dependent ionization leads to 
lower average temperatures inside the flux tubes compared 
to time-independent ionization. However, the average tem- 
peratures inside of flux tubes in models with or without 
time-dependent ionization are nonetheless noticeably higher 
than in the acoustically heated external atmosphere, point- 
ing to the significance of magnetic heating. 

We also pursued corresponding sets of models based on 
LTW frequency spectra instead of monochromatic waves. 
In this type of models the attained time-averaged temper- 
atures were unrealistically low, particularly in models of 
time-dependent ionization. In this case, the combined effects 



of time-dependent hydrogen ionization and of the dynami- 
cal structure related to the propagation of spectral waves 
lead to very strong shocks, which invoked amplified energy 
losses behind the shocks as well as strong atmospheric ex- 
pansions associated with qu a si-adi abati c cooling; see, e.g., 
iGail. Cuntz fc Ulmschneideil (|l990l ) andU oninx fc Pijpers I 
( 1992) for theoretical work on the build-up of wave pressure 
due to shocks. The impact of very strong shocks entailed 
that the time-averaged temperatures are lowest in models 
with small tube spreading, i.e., high magnetic filling factors 
(e.g., case of / = 0.1%) if time-dependent ionization of hy- 
drogen is taken into account. 

Finally, we studied the behaviour of the height- 
dependent time-averaged hydrogen ionization in our set of 
models (see Figs. 8 and 9). We found that the hydrogen ion- 
ization degrees steadily increase at heights beyond 1200 km 
regardless of the magnetic filling factor, i.e., the tube top 
opening radius. However, in general, they remain relatively 
low (i.e., below 40%) in models with time-dependent ion- 
ization (TDI), whereas they are found to be relatively high 
(i.e., between 70% and 90% at heights beyond 2000 km) in 
models without time-dependent ionization (nTDI). This is 
consistent with the general behaviour of time-dependent hy- 
drogen ionization as described, as in this case lower degrees 
of ionization occur at most heights, including regions behind 
shocks, including strong shocks owing to the time-delay of 
ionization in accord to the time-dependent statistical rate 
equations (see Sect. 2.1). 

We also evaluated time-dependent effects for the ion- 
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filling factors of / = 0.1%, 0.01%, and 0.001%. 



ization of calcium and magnesium (see Fig. 8), noting that 
Ca II and Mg II are of significant importance for the facil- 
itation of radiative cooling (see Sect. 2.2). Our simulations 
show that most of the magnesium is ionized to Mg II at most 
heights regardless of the geometrical filling factor of the re- 
spective flux tube. However, for Ca II the height-dependent 
behaviour of the ionization degree is strongly impacted by 
the geometrical spread of the flux tubes (see Fig. 2). Gener- 
ally, for smaller values of /, associated with a lesser dilution 
of the wave energy flux, a larger fraction of calcium is ionized 
to Ca II. The appropriate explanation is the lower dilution 
of the wave energy flux in those models, leading to more 
effective shock heating and higher overall temperatures. 



4 SUMMARY AND CONCLUSIONS 

We pursued sets of time-dependent model simulations for 
solar magnetic flux tubes with focus on longitudinal tube 
waves with and without the consideration of time-dependent 



ionization, particularly pertaining to hydrogen. The treat- 
ment of longitudinal tube waves is motivated by a large va- 
riety o f studies including work by Hasa n fc van BallegooiierJ 
(2008) arguing that longitudinal tube waves are able to pro- 
vide quasi-steady heating sufficient to explain the bright so- 
lar network grains observed in Ca II H and K. We stud- 
ied the dynamics of flux tubes with different magnetic fill- 
ing factors, i.e., 0.1%, 0.01%, and 0.001%, corresponding to 
different tube top opening radii, g uided and motivated by 
previou s obser vational results (e.g.. lSpruit fc Robertsl . ll983l ; 
ISolankil . ll996h . It was found that the dynamic tube struc- 
tures are determined by a complex interplay between shock 
heating, radiative and hydrodynamic cooling (with the lat- 
ter caused by wave pressure), dilution of the wave energy 
flux due to the flux tube geometry, and time-dependent ion- 
ization. We identified pronounced differences between the 
tubes models especially in regard to tubes with different 
tube opening radii and in response to the inclusion or omis- 
sion of time-dependent ionization. 

Similar to the case of acoustic waves, time-dependent 
ionization in LTW models leads to a large range of phe- 
nomena, such as over- and underionization of the flow, in- 
creased temperature jumps at shocks, and modified mean 
(i.e., time-averaged) temperature, density and ionization 
structures. The impact of time-dependent treatment of ion- 
ization for flux tubes becomes evident through comparisons 
of mean temperatures between the different types of models. 
Regarding monochromatic models it was found that time- 
dependent ionization leads to lower average temperatures 
inside flux tubes compared to time-independent ionization, 
an effect already known for acoustically heated models. 

However, the average temperatures inside of flux tubes 
in models with and without time-dependent ionization 
are nonetheless noticeably higher than in the acoustically 
heated external atmosphere. If monochromatic longitudi- 
nal tube waves are used, the mean (i.e., time-averaged) 
temperatures of the tubes generally show a quasi-steady 
increase with height. However, if frequency spectra are 
used, the mean temperatures, in essence, show no rise 
with height at all. This behaviour is initiated by the for- 
mation of very strong shocks (especially in narrow tubes, 
corresponding to small magnetic filling factors) resulting 
in large-scale quasi-adiabatic cooling, leading to unreal- 
istically low temperatures. This type of result has al- 
ready previously been foun d in corresponding 1-D mod- 
els of acoustic waves (e.g., Carlsson fc Steinl Il992l . Il995l ; 
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iRammacher fc Ulmschneider 20031), a beh aviour deemed 
highly unrealistic l|ljlmschneider et all 120051 ). 

An appropriate enhancement of our current models will 
be the calculation of self-consistent 3-D magnetohydrodyn- 
mic (MHD) models. Those models face, however, the prin- 
cipal challenge of the necessity to include both detailed 
multi-level 3-D radiative transfer and 3-D flows including 
the detailed treatment of shock formation and shock in- 
teraction. Importa nt progress has already been made (e.g., 
IStein et al. I, l2009al lbl). but further efforts are needed to ob- 
tain of a concise picture. Tentative insights into the prin- 
cipal properties of these types of future models can be 
attained through inspecting existing 3-D time-dependent 
non-magnetic hydrodynamics models, which also consider 
time -dependent non-equilibrium effects ca used by hydro- 
gen l|Leenaarts fc Wedemever-Bohml [2006), albeit various 
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restrictive assumptions including (but not limited to) the 
lack of back-coupling of the ionization to the equation of 
state. In this type of models it is found that the build- 
up of strong shocks due to shock interaction is largely ab- 
sent, resulting in a lack of unrealistically high cooling be- 
hind the strong shocks pre viously also referred to as "hy- 
drodynamic refrigeration" (|Cuntz fe Muchmorel . 1 1994 ). In 
this case a quasi-steady rise of temperature with height is 
attained, which appears to be in close resemblance to em- 
pirica l solar chromosphere models (e.g.. lAnderson fc Athavl 
1989). Nonetheless, our results based on time-dependent ion- 
ization and LTW frequency spectra allow insight into the 
limiting case of 1-D geometry, while also noting that our 
models based on time-dependent ionization and monochro- 
matic waves are expected to be approximately reflective of 
physical reality. 

An alternative, or perhaps supplementary, way of sup- 
plying chromospheric heating might be given through am- 
bipola r diffusion as described by iKhomenko fc Colladosl 
l|2012h . Here the presence of neutrals, together with the de- 
crease with height of the collisional coupling, leads to devi- 
ations from the classical ma gnetohydrodynamic behaviou r 
of the chromospheric plasma. iKhomenko fc Colladosl (|2012l ) 
pointed out that a relative net motion occurs between the 
neutral and ionized components, referred to as ambipolar 
diffusion. According to this model, the dissipation of cur- 
rents in the chromosphere is enhanced by orders of magni- 
tude due to the action of ambipolar diffusion, as compared 
with the standard ohmic diffusion. The authors proposed 
that a significant amount of magnetic energy can be released 
to the chromosphere just by existing force-free 10-40 G mag- 
netic fields there. 

Additional stud ies were given by 

iFedun. Erdelyi fc Shelvasl (120091'). I Vigeesh. Hasan fc Steinerl 
l|2009h . and lErdelvi fc Fedunl (|201ol ). iFedun et all (<2009r ) 
studied the oscillatory response of the 3-D solar photo- 
sphere to the leakage of photospheric motion. They found, 
among other results, that high-frequency waves propagate 
from the lower atmosphere across the transition region 
experiencing relatively low reflection, and transmitting 
most of their energy into the corona, and, furthermore, 
that the magnetic field acts as a waveguide for both high- 
and low-frequency waves originating from th e photosphere 
and p ropagating up into the solar corona. IVigeesh et al.l 
(2009) provided a targeted study on wave propagation 
and energy transport in the magnetic network of the 
Sun based on 2-D MHD simulations, which among other 
results identified the l imited capacity of acoustic waves. 
lErdelvi fc Fedunl (|2010h investigated the oscillatory modes 
of a magnetically twisted compressible flux tube embedded 
in a compressible magnetic setting, including applications 
to solar magneto-seismology. 

As an overarching statement concerning our study we 
conclude that the significance of time-dependent ionization 
identified in the simulations of longitudinal tube waves is a 
stark motivation to also consider this type of effect in fu- 
ture models of transverse and torsional tube waves. This 
will allow to obtain a more detailed picture of the dynamics 
and energetics of solar- type outer atmospheres. This obvi- 
ous suggestion is also supported by the repeatedly obtained 
finding that longitudinal flux tube waves, as gauged through 
models considering time-dependent ionization phenomena as 



done in the present study, are insufncent to supply an ade- 
qute amount of energy for balancing cor onal he a ting i n the 
view of early and updated estimates by iGudell l|2007h and 
others. 
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